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ABSTRACT 

This  paper  derives  a  system  of  equations  that  approximately  govern  small-amplitude  perturbations  in  a 
nonprecipitating  cloudy  vortex.  The  cloud  coverage  can  be  partial  or  complete.  The  model  is  used  to 
examine  moist  vortex  Rossby  wave  dynamics  analytically  and  computationally.  One  example  shows  that 
clouds  can  slow  the  growth  of  phase-locked  counter-propagating  vortex  Rossby  waves  in  the  eyewall  of  a 
hurricane-like  vortex.  Another  example  shows  that  clouds  can  (indirectly)  damp  discrete  vortex  Rossby 
waves  that  would  otherwise  grow  and  excite  spiral  inertia-gravity  wave  radiation  from  a  monotonic  cyclone 
at  high  Rossby  number. 


1.  Introduction 

There  have  been  many  attempts  to  understand  how 
wave  motions  within  a  tropical  cyclone  regulate  its  in¬ 
tensity  (see  below).  Most  of  these  efforts  were  based  on 
wave  theories  that  neglect  moisture.  Since  clouds  per¬ 
vade  tropical  cyclones,  such  an  approximation  is  not 
soundly  justified.  The  warming  or  cooling  of  air  parcels 
by  condensation  or  evaporation  of  water  content  can 
dramatically  affect  the  wave  dynamics.  In  this  paper,  we 
will  begin  to  see  how. 

In  particular,  we  will  examine  waves  in  a  nonprecipi¬ 
tating  cloudy  vortex  whose  basic  state  has  no  secondary 
circulation  and  is  gravitationally  stable.  In  general,  the 
basic  state  has  fluctuations  of  cloudy  and  unsaturated 
air,  but  those  fluctuations  are  assumed  to  occur  on 
smaller  scales  than  the  waves  of  interest.  We  will  not 
consider  the  effects  of  boundary  layer  fluxes,  radiative 
heat  fluxes,  nor  eddy  diffusion  on  the  wave  dynamics. 
Furthermore,  we  will  suppose  that  perturbations  are 
hydrostatic  and  governed  by  linearized  equations. 

In  sections  2-5  we  will  formally  develop  this  model. 
Although  it  applies  to  both  hydrostatic  inertia-gravity 
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(IG)  waves  and  vortex  Rossby  (VR)  waves,  we  will 
focus  our  attention  on  the  latter. 

In  section  6,  we  will  briefly  reexamine  the  eyewall 
instability  of  a  baroclinic  cyclone  that  resembles  a  cat¬ 
egory  3  hurricane.  Nolan  and  Montgomery  (2002)  re¬ 
ported  that  this  particular  vortex  destabilizes  most  rap¬ 
idly  by  the  interaction  of  counter  propagating  VR 
waves  on  opposite  edges  of  the  eyewall.  We  will  show 
that  increased  fractional  cloud  coverage  slows  the  in¬ 
stability,  and  thereby  hinders  asymmetric  eyewall 
breakdown.  The  consequence  on  hurricane  intensity  is 
a  topic  of  ongoing  study.  On  the  one  hand,  potential 
vorticity  mixing  after  eyewall  breakdown  can  directly 
reduce  the  maximum  tangential  wind  speed  of  a  hurri¬ 
cane  (e.g.,  Schubert  et  al.  1999;  Kossin  and  Schubert 
2001).  On  the  other  hand,  the  thermo-fluid  dynamics 
connected  to  eyewall  breakdown  is  subtle,  and  may  ac¬ 
tually  lead  to  a  stronger  hurricane  (Emanuel  1997; 
Montgomery  et  al.  2002;  Persing  and  Montgomery 
2003;  Montgomery  et  al.  2006). 

In  section  8,  we  will  examine  how  moisture  influ¬ 
ences  the  discrete  VR  waves  of  a  barotropic  cyclone, 
whose  vertical  vorticity  decreases  monotonically  with 
radius.  These  waves  account  for  precessing  tilts  and 
elliptical  (triangular,  square,  etc.)  deformations  of  the 
vortex  core.  If  the  Rossby  number  of  the  cyclone  ex¬ 
ceeds  unity,  its  baroclinic  VR  waves  can  efficiently  ex- 
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cite  frequency-matched  spiral  IG  waves  in  the  environ¬ 
ment  (Ford  1994a, b;  Polvani  et  al.  1994;  Plougonven 
and  Zeitlin  2002;  Schecter  and  Montgomery  2003,  2004, 
2006).  Surprisingly,  producing  IG  wave  radiation  com¬ 
pels  a  VR  wave  to  grow.  Such  growth  is  usually  op¬ 
posed  by  a  critical  layer  that  resonantly  absorbs  wave 
activity  (Briggs  et  al.  1970;  Pillai  and  Gould  1994; 
Schecter  et  al.  2000,  2002;  Balmforth  et  al.  2001;  Mallen 
et  al.  2005).  We  will  show  that  moisture  tends  to  en¬ 
hance  damping  by  the  critical  layer  and  thereby  sup¬ 
press  the  radiative  instability.  Consequently,  cloud  cov¬ 
erage  can  inhibit  the  loss  of  angular  momentum  by  IG 
wave  radiation.  Such  loss  might  otherwise  merit  some 
consideration  in  the  angular  momentum  budget  of  a 
hurricane  (e.g.,  Chow  and  Chan  2003). 

The  theory  of  discrete  VR  waves,  especially  of  the  tilt 
mode,  also  provides  a  framework  for  understanding  the 
response  of  a  vortex  to  ambient  vertical  shear.  Recent 
studies  (e.g.,  Jones  1995;  Reasor  and  Montgomery 
2001;  Schecter  et  al.  2002;  Reasor  et  al.  2004)  have 
shown  that  reducing  the  static  stability  of  the  atmo¬ 
sphere  typically  increases  (i)  the  azimuthal  phase  ve¬ 
locity  and  (ii)  the  decay  rate  of  the  discrete  tilt  mode  of 
a  monotonic  (or  quasi-monotonic)  cyclone.  Both  of 
these  modifications  increase  the  cyclone’s  natural  resis¬ 
tance  to  misalignment  under  the  influence  of  ambient 
shear.  As  we  will  show  in  section  8,  increasing  the  cloud 
coverage  has  similar  effects  on  the  tilt  mode  and 
thereby  also  improves  vortex  resilience  (cf.  Patra  2004). 
This  result  is  quite  reasonable,  in  light  of  the  aforemen¬ 
tioned  studies,  if  one  connects  denser  cloud  coverage  to 
lower  static  stability  (e.g.,  Durran  and  Klemp  1982). 

Another  class  of  VR  waves  comprises  continuum 
mode  disturbances  that  are  sheared  by  differential  ro¬ 
tation  of  the  mean  flow.  Sheared  VR  waves  that  are 
generated  by  localized  bursts  of  convection  can  trans¬ 
port  angular  momentum  toward  the  radius  of  maximum 
wind,  and  thereby  intensify  a  cyclone  (e.g.,  Montgom¬ 
ery  and  Kallenbach  1997;  Montgomery  and  Enagonio 
1998;  Moller  and  Montgomery  1999,  2000;  Enagonio 
and  Montgomery  2001).  They  may  also  have  some  con¬ 
nection  to  moving  spiral  rainbands  (e.g.,  Montgomery 
and  Kallenbach  1997).  Although  we  will  not  directly 
examine  the  influence  of  clouds  on  sheared  VR  waves, 
a  theory  based  on  our  moist  model  would  be  a  modest 
formal  extension  of  the  existing  dry  theory. 

On  a  technical  note,  the  moist  wave  equations  of 
section  4  are  nearly  identical  to  the  linearized  dry  hy¬ 
drostatic  primitive  equations.  Moisture  simply  reduces 
the  Eulerian  rate  of  cooling/warming  by  updrafts/down¬ 
drafts  by  a  factor  Y  [see  Eq.  (31)].  The  Y  distribution  is 
a  property  of  the  basic  state  of  the  vortex,  and  has  only 
radial  and  vertical  variation.  Its  value  indicates  the  av¬ 


erage  reduction  of  the  vertical  buoyancy  restoring  force 
of  air  parcels  along  an  azimuthal  circuit  (see  appendix 
G).  As  the  fraction  of  cloudy  air  in  the  circuit  increases 
from  zero,  Y  decreases  from  unity.  In  section  5,  we  will 
present  a  formula  for  Y  [Eq.  (37)]  that  makes  a  baro- 
clinic  vortex  neutral  to  slantwise  convection.  In  addi¬ 
tion,  we  will  examine  the  spatial  distribution  of  Y  for  a 
realistically  simulated  hurricane.  As  expected,  we  will 
find  that  Y  is  smallest  in  the  very  moist  eyewall. 

The  moist  wave  equations  presented  here  are  an  im¬ 
provement  over  the  usual  dry  model,  and  their  solu¬ 
tions  provide  insight  into  the  most  basic  effects  of  cloud 
coverage.  Nevertheless,  they  exclude  several  factors 
(see  paragraph  2)  that  might  significantly  alter  wave 
motions  in  actual  or  realistically  simulated  hurricanes 
(cf.  Wang  2002a,b;  Chen  and  Yau  2001;  Chow  et  al. 
2002;  Chen  et  al.  2003;  etc.).  A  more  comprehensive 
theory  should  be  the  goal  of  future  research. 

The  remainder  of  this  paper  is  organized  as  follows: 
section  2  describes  our  parameterization  of  moist  pro¬ 
cesses  in  a  cloudy  vortex.  Section  3  presents  the  non¬ 
linear  hydrostatic  fluid  equations  that  are  supposed  to 
govern  the  vortex.  Section  4  presents  the  linearized 
wave  equations.  Section  5  discusses  some  important  Y 
distributions.  Section  6  illustrates  how  cloud  coverage 
can  temper  an  asymmetric  eyewall  instability.  Section  7 
contains  a  formal  discussion  of  the  eigenmode  problem 
for  a  cloudy  vortex.  Section  8  examines  how  cloud  cov¬ 
erage  can  cause  or  accelerate  the  decay  of  discrete  VR 
waves  in  a  monotonic  cyclone.  Such  decay  inhibits  the 
radiation  of  spiral  IG  waves  and  improves  vortex  resil¬ 
ience.  Section  9  summarizes  our  main  results.  Appen¬ 
dices  A-F  provide  some  analytical  and  computational 
details  for  those  who  may  wish  to  continue  the  analysis 
of  waves  in  a  cloudy  vortex  beyond  the  scope  of  this 
paper.  Appendix  G  informally  discusses  the  physical 
meaning  of  Y. 

2.  Moisture  parameterization 

a.  The  water  mass  equations 

Suppose  that  water  has  only  two  phases:  vapor  (v) 
and  liquid  (/).  Further  suppose  that  there  is  no  precipi¬ 
tation,  and  that  there  is  no  diffusion  of  water  mass. 
Then,  the  total  water  mixing  ratio, 

qt  =  q^  +  qo  (i) 

is  conserved  along  material  trajectories.  That  is. 


The  water  vapor  mixing  ratio  is  here  determined  by 
the  diagnostic  equation, 
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in  which  p)  is  the  saturation  vapor  mixing  ratio, 

T  is  temperature,  and  p  is  total  pressure. 

b.  The  heat  equation 

The  entropy  of  a  moist  fluid  parcel,  per  unit  mass  of 
dry  air,  is  given  by  [e.g.,  Emanuel  1994,  Eq.  (4.5.9)] 


■5m  =  (Cpd  +  q,Ci)lniT)  -  RdMp^) 

+  KqJT  -  q^RM(e/e^). 


We  will  assume  that 


dSm 

dt 


=  0. 


(4) 


(5) 


P  ( Po\^‘^^‘^P‘‘ 


(6) 


0p=  ^ 


(7) 


h“p{p,s^,q,)  q,<q^.^, 

epip,s^,q,)  q,>q^^. 

For  convenience,  we  introduce  the  saturation  variable 

~  qt  qvii:^  (^) 

which  is  positive  (negative)  for  saturated  (unsaturated) 
parcels.  In  addition,  let  //(p)  denote  the  unit  step  func¬ 
tion.  Then,  we  may  rewrite  Eq.  (7)  as  follows: 


2(X)- 


275 


285  295 

T(K) 


— 1(1- 
.405 


Above,  Cpj  and  C/  are  the  specific  heats  of  dry  air  (at 
constant  pressure)  and  of  liquid  water;  R^  and  R^  are 
the  gas  constants  of  dry  air  and  of  water  vapor;  p^  is  the 
pressure  of  dry  air,  whereas  e  and  e*(T)  are  the  actual 
and  saturation  values  of  the  water  vapor  pressure; 
LJJ")  is  the  latent  heat  of  vaporization. 

We  may  now  derive  an  equation  for  the  density  po¬ 
tential  temperature,  which  is  defined  by 


Here,  p  is  the  total  mass  density  of  a  fluid  parcel  and  p^ 
is  a  reference  pressure,  taken  to  be  10®  Pa. 

In  the  present  model,  three  variables  specify  the  ther¬ 
modynamic  state  of  a  moist  air  parcel  [Eq.  (3)  elimi¬ 
nates  gj.  We  may  choose  these  variables  to  be  p,  s^, 
and  q,.  In  terms  of  these  variables,  the  functional  form 
of  Op  differs  for  saturated  (s)  and  unsaturated  (u)  air 
parcels;  that  is. 


Fig.  1.  Contour  plot  of  minus  x  for  p,  =  0.015.  The  contour 
labels  are  in  units  of  K  Pa“^.  Note  that  x  is  discontinuous  across 
the  curve  p)  =  <?r,  which  separates  (left)  saturated  and 

(right)  unsaturated  regions.  The  dashed  curves  are  contours  of 
constant  moist  entropy  s^.  The  value  of  entropy  increases  with 
the  superscript  i. 

The  material  derivative  of  density  potential  tempera¬ 
ture  has  three  terms;  specifically, 

c^^dp/^\  ^  dqJdOp 

dt  ~  dt  dt  \ds„,)p^^  dt  \dq,,,^^p 

(10) 

Inserting  Eqs.  (2),  (5),  and  (9)  into  the  right-hand  side 
of  Eq.  (10)  yields 


dOp  dp 
dt  dt 


(  Sq^, 


in  which 


,  dOp 


(11) 


(12) 


Here,  we  have  used  the  fact  that  dHIdp.  =  5(p.),  in 
which  S(p,)  is  the  Dirac  distribution.  Figure  1  is  a  con¬ 
tour  plot  of  minus  x  for  an  air  parcel  with  q,  =  0.015. 
This  value  of  q^  is  characteristic  of,  say,  a  parcel  in  the 
lower  region  of  a  tropical  storm.  Note  that  the  negative 
value  of  X  increases  by  orders  of  magnitude  as  a  parcel 
moves  from  an  unsaturated  state  (p,  <  0)  to  a  saturated 
state  (p  >  0)  along  any  path,  such  as  a  moist  adiabat 
(dashed  curve). 

We  may  simplify  Eq.  (11)  by  observing  that 


(9) 


6“  =  01  if  p  =  0. 


Op  =  //(-p)0“  +  H{p.)(fp. 


(13) 
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Therefore,  Eq.  (11)  reduces  to 


ddp  dp 

m  ^  ^~dt  ^ 


(14) 


which  constitutes  our  working  heat  equation.  In  a  dry 
atmosphere,  x  =  0  6p  is  conserved. 

In  appendix  A,  we  derive  the  following  equations  for 


the  thermodynamic  derivatives  that  appear  in  the  defi¬ 
nition  of  X- 


dp 


TRd 

PCpd 

X  I  1  - 


\p  )  1  +  ?( 


1  +  qi/e 
1  +  qtCpJCpd 


(15) 


and 


\  ^P  )s„,q,  PCpd  \  P  /  1  + 


1  +  ■ 


R.e 


(1  +  1  + 


RdT 


4  ^v:i:^pv^^pd  P^l^^pd 


R^Cp^T 


■  (16) 


Here,  we  have  introduced  the  ratio  of  gas  constants  e  = 
RJRv  =  0.622.  Since  =  q^^{T,  p)  and  =  LJT), 
the  right-hand  sides  of  Eqs.  (15)  and  (16)  are  functions 
of  the  measurable  quantities  T,  p,  and  q,. 


d(t>  Rd  (  P 

Sp  P  Uo/ 

and  the  mass  continuity  equation  yields 


(20) 


3.  The  moist  primitive  equations 

For  the  remainder  of  this  paper,  we  will  assume  that 
vertical  accelerations  are  small  and  that  hydrostatic  bal¬ 
ance  applies.  Consequently,  we  may  use  p  as  the  verti¬ 
cal  coordinate  instead  of  height.  In  addition,  we  will 
assume  that  the  vortex  occupies  a  sufficiently  small  sec¬ 
tion  of  the  globe  where  the  /-plane  approximation  is 
valid. 

In  pressure  coordinates,  the  horizontal  velocity  equa¬ 
tion  takes  the  form 

Id  \  au 

+  +  +  (17) 

in  which  u  is  the  horizontal  velocity  vector,  w  =  dpidt, 
f  is  the  constant  Coriolis  parameter,  and  (p  =  gz  is  the 
geopotential.  Furthermore,  z  is  the  vertical  unit  vector, 
and  Vp  is  the  horizontal  gradient  operator,  with  spatial 
derivatives  evaluated  at  constant  p. 

The  heat  equation  (14)  and  conservation  of  water 
mass  (2)  take  the  forms 

+  »  (18) 

and 

/  d  \  dq^ 

[-  +  U-Vp)q,=  -a.-,  (19) 

respectively. 

In  addition  to  the  above  prognostic  equations,  there 
are  two  diagnostic  equations.  Hydrostatic  balance 
yields 


dw 

dp 


=  -v„ 


u. 


(21) 


The  auxiliary  equation  for  hydrostatic  potential  vor- 
ticity. 


n  -  r,  •  V0p,  (22) 

is  given  by 

/  d  \  an 

+  +  =  (23) 

in  which  V  =  —  zd/dp,  and 

K)  =  V  X  u  -f  /z.  (24) 

Note  that  we  have  used  a  convention  in  which  the  hori¬ 
zontal  and  vertical  components  of  V  have  different  di¬ 
mensions  (inverse  length  and  inverse  pressure).  We 
have  used  a  similar  convention  for  tj.  Of  course  each 
term  of  tj  •  V  (for  example)  has  the  same  dimension. 


4.  The  wave  equations 

This  section  presents  a  linear  model  for  wave  dynam¬ 
ics  in  a  cloudy  vortex  whose  basic  state  has  no  second¬ 
ary  circulation.  We  will  use  a  vortex-centered  cylindri¬ 
cal  coordinate  system  in  which  r  and  <p  are  the  radius 
and  azimuth,  and  u  and  v  are  the  radial  and  azimuthal 
velocities.  Overbars  and  primes  will  mark  basic  state 
and  perturbation  fields,  respectively. 
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a.  Balanced  basic  state 


Furthermore,  the  heat  equation  (18)  may  be  written 


The  basic  state  of  an  atmospheric  vortex  is  charac¬ 
terized  by  its  azimuthal  velocity  field  v{r,  p).  We  may 
define  u  to  be  a  tempero-azimuthal  mean  of  v.  For  no- 
tational  convenience,  we  introduce  the  following  auxil¬ 
iary  fields: 

VL{r,p)  =  v/r,  ^{r,  p)  =  r~^d{rv}/dr 

V(r,p)  =  l  +  f,  I(C  p)  =  2CL+f.  ^  ^ 

Here,  fl  is  the  angular  rotation  frequency,  |  is  the  rela¬ 
tive  vertical  vorticity,  f)  is  the  absolute  vertical  vorticity, 
and  I  is  the  modified  Coriolis  parameter.  By  assump¬ 
tion,  the  basic  state  has  no  secondary  circulation,  that  is, 
cd,  M  =  0. 

The  basic-state  geopotential  p)  and  density  po¬ 
tential  temperature  Qjj,  p)  are  defined  to  satisfy  gra¬ 
dient  balance. 


dr 


(26) 


and  hydrostatic  balance. 


d(j)  Rd  I  P 

dp  P  \po) 


(27) 


Taking  the  partial  derivative  of  Eq.  (26)  with  respect  to 
pressure,  making  use  of  Eq.  (27),  and  simplifying,  we 
obtain  the  thermal  wind  relation. 


ddp  P  ( dv 

dr  Rd\p  /  ^  dp' 


(28) 


The  basic  state  also  has  a  water  mass  profile  q,{r,  p), 
which  we  take  to  be  a  tempero-azimuthal  mean  of  q,. 
We  further  define  )((r,  p)  to  be  a  tempero-azimuthal 
mean  of  x-  Note  that  the  difference  between  )((r,  p)  and 
x(0p,  q,,  p)  may  be  substantial  if  q,  fluctuates  between 
saturated  and  unsaturated  values  within  the  t  and  (p 
averaging  intervals. 


b.  Small  perturbations 

Each  field  that  describes  the  vortex  is,  of  course,  the 
sum  of  its  basic  state  and  a  time-dependent  perturba¬ 
tion.  Neglecting  terms  that  are  quadratic  in  the  pertur¬ 
bation  fields,  the  two  components  of  the  horizontal  ve¬ 
locity  equation  (17)  may  be  written 


and 


u'  = 


d^' 

dr 


(29) 


d  _  d  _ 

—  -f  fi  —  In'  =  —r\u' 


dt 


dip 


dv  1  dip' 

dp  r  dip  ' 


(30) 


dOp 


(31) 


in  which  we  have  introduced  the  buoyancy  reduction 
factor. 


r(r,p) 


dOp/dp  -  X 
dOp/dp 


(32) 


Section  5  discusses  Y  distributions  that  might  have  spe¬ 
cial  relevance  to  atmospheric  vortices.  Appendix  G  ex¬ 
plains  in  detail  how  Y  measures  a  local  reduction  of  the 
vertical  buoyancy  restoring  force  that  acts  on  vertically 
displaced  air  parcels. 

In  addition  to  the  above  prognostic  equations,  we 
have,  from  hydrostatic  balance. 


dp~  P  \po) 
and,  from  mass  continuity, 

db)'  1  d{ru')  1  dv' 

dp  r  dr  r  dip  ' 


(33) 


(34) 


We  need  not  consider  the  linearized  equation  for  total 
water  mass,  since  the  above  system  (29)-(34)  does  not 
depend  explicitly  on  q',,  and  is  therefore  closed. 

Note  that  the  above  wave  equations  are  the  same  as 
those  for  a  dry  vortex,  with  dOpIdp  — >  Yddpidp. 

Note  also  that  linearization  of  the  heat  equation  is  a 
bold  approximation.  Linear  theory  assumes  that  qua¬ 
dratic  perturbation  terms  are  subdominant  to  linear 
terms  in  the  equations  of  motion.  However,  co'x'  is  not 
necessarily  negligible.  Since  x  varies  discontinuously 
between  unsaturated  and  cloudy  air,  x’  can  have  fluc¬ 
tuations  of  order  x-  If  these  fluctuations  have  small 
spatiotemporal  scales,  they  may  not  appreciably  influ¬ 
ence  the  large-scale  waves  of  interest.  We  defer  testing 
the  accuracy  of  linearization,  with  high-resolution  non¬ 
linear  simulations,  until  a  future  time. 


5.  The  buoyancy  reduction  factor 

For  the  remainder  of  this  paper,  we  will  treat  the 
buoyancy  reduction  factor  Y(r,  p)  as  an  adjustable  pa¬ 
rameter.  We  will  attempt  to  choose  Y  distributions  that 
have  some  meteorological  significance. 

a.  Y  for  symmetric  stability 

There  is  a  common  belief  that  intense  atmospheric 
vortices,  hurricanes  in  particular,  have  basic  states  that 
are  nearly  neutral  to  slantwise  symmetric  perturbations 
(e.g.,  Emanuel  1986).  A  necessary  condition  for  the 
symmetric  stability  of  a  cloudy  baroclinic  vortex,  with 
no  secondary  circulation,  may  be  written  as  a  local  con¬ 
straint  on  Y. 
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Fig.  2.  The  (left)  azimuthal  velocity,  (middle)  vertical  velocity,  and  (right)  buoyancy  reduction  factor  of  a  mature  hurricane,  simulated 
with  RAMS.  The  dashed  curves  in  the  vertical  velocity  plot  are  zero  contours.  In  all  plots,  the  radius  and  height  axes  are  scaled  linearly. 
The  data  set  used  to  produce  this  figure  was  kindly  provided  by  Dr.  M.  Nicholls  and  Dr.  J.  Persing. 


To  begin  with,  let 


d  1  d(ru')  du' 

—  ^0  and - ^  — 

dtp  r  dr  dr 

in  Eqs.  (29)-(31)  and  (34).  In  addition,  let 

.  _  dv  dO  dd 

I,  T],  —  ,  — —  ,  Y  — ^  constants. 


(35) 


(36) 


along  with  the  coefficient  of  Op  on  the  right-hand  side  of 
Eq.  (33).  The  above  substitutions  are  justifiable  only  if 
(i)  the  perturbation  is  roughly  axisymmetric,  (ii)  the 
perturbation  fields  vary  much  more  rapidly  than  the 
mean  fields  with  r  and  p,  and  (iii)  the  radial  wavelength 
of  the  perturbation  is  much  less  than  r.  We  will  refer  to 
the  simplified  equations  as  the  local  axisymmetric 
model. 

Appendix  B  shows  that  perturbations  are  stable  in 
the  local  axisymmetric  model  if  and  only  if 


Y>Y„,(r,/7)- 
or  equivalently. 


[^{dv/dp)f  /  ^  \ 

Rapad0p/dp)^[p  ) 


(38) 


in  which  the  pressure  derivative  is  along  a  contour  of 
constant  angular  momentum,  M  =  rv  +  fr^/2.  For  glob¬ 
al  stability,  this  condition  must  hold  at  every  point  in 
the  r-p  plane  (cf.  Ooyama  1966). 

We  have  not  proven  that  condition  (37)  guarantees 
nongrowth  of  vortex-scale  symmetric  perturbations. 
Nevertheless,  if  Y  =  Y„„  we  will  casually  say  that  the 
vortex  is  neutral  to  slantwise  convection.  In  section  6, 
we  will  see  how  the  growth  of  an  asymmetric  perturba¬ 
tion  can  change  as  a  vortex  moistens  and  becomes  neu¬ 
tral  to  slantwise  convection. 


It  is  worth  mentioning  that  even  complete  cloud  cov¬ 
erage  can  leave  sections  of  a  vortex  symmetrically 
stable.  Therefore,  although  setting  Y  =  Y„,  is  of  theo¬ 
retical  interest,  it  is  not  always  physical. 

b.  Y  for  a  simulated  hurricane 

Let  us  now  consider  the  Y  distribution  of  a  mature 
hurricane  in  a  “cloud-resolving”  numerical  simulation. 
This  particular  simulation  was  performed  with  the  Re¬ 
gional  Atmospheric  Modeling  System  (RAMS),  devel¬ 
oped  at  Colorado  State  University  (see  Cotton  et  al. 
2003).  It  was  optionally  set  to  include  liquid  cloud  and 
rain,  but  not  ice. 

Figure  2  shows  the  instantaneous  azimuthal  means  of 
the  tangential  and  vertical  velocity  fields  of  the  simu¬ 
lated  hurricane.  In  addition,  it  shows  the  buoyancy  re¬ 
duction  factor  Y.  To  evaluate  Y  [the  right-hand  side  of 
Eq.  (32)],  we  used  the  instantaneous  azimuthal  means 
of  Op  and  To  evaluate  x  at  each  grid  point  we  used 
Eq.  (15),  Eq.  (16),  and  the  RAMS  approximation  for 
q^-AT,  p)-  As  expected,  Y  is  minimal  in  the  main  up¬ 
draft  region,  where  the  air  is  very  moist.  This  region 
also  contains  the  strongest  tangential  winds.  Low  Y  at 
high  u  is  a  key  result  to  extract  from  this  figure. 

Although  tempting,  we  cannot  use  the  theory  of  sec¬ 
tion  5a  to  properly  evaluate  the  symmetric  stability  of 
the  RAMS  hurricane.  Condition  (37)  strictly  applies  to 
a  vortex  that  has  no  mean  secondary  circulation.  In 
contrast,  the  RAMS  hurricane  is  maintained  by  a 
strong  radial  inflow  (not  shown)  that  has  a  peak  veloc¬ 
ity  of  about  30  m  s"’^  near  the  surface.  Still,  the  negative 


^  Because  of  strong  secondary  circulation,  the  azimuthal  mean 
of  0p,  averaged  over  any  time  interval,  does  not  exactly  satisfy  the 
thermal  wind  equation  (28).  Therefore,  the  plotted  buoyancy  re¬ 
duction  factor  differs  slightly  from  the  formally  defined  T. 
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values  of  Y  at  the  base  of  the  main  updraft  might  indi¬ 
cate  that  the  vortex  is  gravitationally  unstable  in  that 
region. 

6.  Eyewall  instability 

The  circular  eyewall  of  a  mature  hurricane  can  de¬ 
stabilize  and  evolve  into  a  ring  of  mesovortices,  which 
itself  can  destabilize  and  rapidly  spread  angular  mo¬ 
mentum  into  the  eye.  This  process,  the  asymmetric  in¬ 
stability  of  a  cylindrical  vortex  sheet,  has  been  studied 
extensively  in  laboratory  experiments  (e.g.,  Weske  and 
Rankin  1963;  Ward  1972;  Vladimirov  and  Tarasov 
1980;  Montgomery  et  al.  2002).  It  has  also  been  studied 
theoretically/computationally  without  much  consider¬ 
ation  of  moisture  (e.g.,  Michalke  and  Timme  1967;  Ro- 
tunno  1978;  Gall  1983;  Davidson  1990,  chap.  6;  Schu¬ 
bert  et  al.  1999;  Kossin  and  Schubert  2001;  Nolan  and 
Montgomery  2002;  Terwey  and  Montgomery  2002). 
However,  hurricane  eyewalls  are  covered  with  rela¬ 
tively  deep  clouds  (see  Fig.  2).  In  this  section,  we  will 
examine  a  case  in  which  such  cloud  coverage  gravely 
decelerates  the  linear  stage  of  an  asymmetric  eyewall 
instability. 


coupled  to  an  invertibility  relation, 


h  d  d  I  l\ 


(41) 


The  length  scale  Ij^  that  appears  in  the  differential  op¬ 
erator  is  the  Rossby  deformation  radius.  It  is  here  de¬ 
fined  by  Ir  -  Vghjf,  in  which  /i„  is  the  ambient  sur¬ 
face  height. 

The  unperturbed  PV  distributions  of  the  rings  under 
consideration  have  the  form 


n,g(r)  =  Zo  +  Z^H{r  -  r,)  -  Z^H{r  -  r^),  (42) 

in  which  <  r2,  Zo  and  Zi  are  positive,  and  Z2  = 
Zq  +  Z^.  As  r  increases,  the  PV  jumps  from  Zq  to  Zq  + 
Zj  at  the  radius  and  then  drops  to  zero  at  rj. 

The  unperturbed  angular  rotation  frequency  is  for¬ 
mally  given  by 

(43) 

The  Green  function  that  appears  in  the  radial  integral  is 
the  n  =  0  member  of  the  family  that  is  defined  by 


a.  Preliminary  consideration:  The  quasigeostrophic 
shallow -water  potential  vorticity  ring 


The  shallow-water  potential  vorticity  ring  is  perhaps 
the  simplest  analogue  of  a  mature  hurricane.  In  this 
analogy,  one  relates  cloud  coverage  to  a  reduction  of 
the  gravitational  restoring  force.  As  a  preparatory  ex¬ 
ercise,  we  here  examine  the  effect  of  lowering  g  (the 
gravitational  acceleration)  on  the  asymmetric  break¬ 
down  of  a  shallow-water  potential  vorticity  ring.  For 
analytical  simplicity,  we  confine  our  discussion  to  the 
quasigeostrophic  (QG)  parameter  regime.  Although 
hurricanes  have  basic  states  that  are  ageostrophic,  we 
may  expect  their  vortex-scale  shear-flow  instabilities  to 
have  QG-like  properties  (e.g.,  Montgomery  and  Sha¬ 
piro  1995).  We  will  return  to  this  matter  in  section  6c. 

In  QG  theory,  fluid  motions  are  assumed  to  be  slow 
compared  to  inertial  oscillations.  Consequently,  the 
dominant  component  of  the  velocity  field  is  the  diver¬ 
gence-free  geostrophic  wind. 


5i/;„  di/;„ 

- +  — s  A 

r  d(p  dr 


(39) 


Here,  we  have  introduced  the  geostrophic  streamfunc- 
tion,  if)g  =  gh'lf,  in  which  h'  is  the  surface-height 
anomaly  of  the  shallow-water  layer. 

The  dynamics  of  a  QG  shallow-water  ring  is  gov¬ 
erned  by  material  conservation  of  PV, 


(40) 


(id  d  n^  1  \ 


-  f) 

r 


(44) 


Requiring  regularity  at  the  origin  and  evanescence  as 
r  ^  00,  we  have 


G„(r|r)  =  -I„{r^/lf,)K„{r^/lj^),  (45) 

in  which  /„  and  are  modified  Bessel  functions,  and 
r>(r<)  is  the  greater  (lesser)  of  r  and  f.  We  note  that  if 
Zq  is  zero,  fi  is  negative  in  the  vortex  core.  Despite  its 
anticyclonic  interior,  this  shear  flow  still  resembles  a 
mature  hurricane  in  that  fi  rapidly  increases  with  radius 
(to  a  positive  value)  between  and  ^2,  and  then  decays 
as  r  tends  to  infinity. 

In  what  follows,  we  will  use  the  notation  flj  and  fl2  to 
denote  the  values  of  fi  at  r^  and  r2,  respectively.  Using 
various  analytical  properties  of  modified  Bessel  func¬ 
tions  that  are,  for  example,  listed  in  Abramowitz  and 
Stegun  (1972),  one  finds  from  Eq.  (43)  that 


ai  =  ZiGi(ri|ri)-Z2-Gi(rik2), 


(46) 


and 


02  =  Zi  -  G,(r,  I  r2)  -  Z2Gi(r2 1  /'2).  (47) 

'2 

A  QG  shallow-water  ring  supports  two  discrete 
modes  for  each  azimuthal  wavenumber  n.  Both  have 
the  form 


dt 


r  dr  d(p 


r  d(p  dr 
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n;^(r,  cp,  t)  =  [niS(r  -  r^)  +  n2S(r  -  +  c.c., 

(48) 

in  which  Hi  and  112  sr®  complex  constants,  and  c.c. 
denotes  the  complex  conjugate  of  the  preceding  term. 
Evidently,  both  modes  contain  a  pair  of  frequency- 
matched  VR  waves  on  opposite  edges  of  the  ring.  The 
eigenfrequency  is  generally  complex;  that  is,  v  =  vj^  + 
ivj,  in  which  Vf,  and  vj  are  real. 

Standard  analysis  leads  to  the  following  modal  dis¬ 
persion  relation: 

v^  =  --±-\/b^-4c,  (49) 


in  which 

b  =  n[ZiG„(G|ri)  -  Z2G„(r2k2)  -  (fli  -f  fia)], 

(50) 

and 

c  =  n\n^fl2  +  Cl^Z2G„{r2\r2)  -  Cl2ZiG„{rj^\rj) 

+  ZiZ2[G„(ri|r2)^  -  G„(rilo)G„(r2|r2)]}.  (51) 

If  «  =  1,  it  can  be  shown  that  c  =  0  and  the  discrete 
modes  are  neutral.  For  «  >  2,  there  will  exist  a  growing 
mode  and  a  conjugate  damped  mode  if  the  discriminant 
b^  —  4c  is  negative. 

With  the  aid  of  symbolic  computation,^  we  have 
found  that 


c 


l%n^{n^  -  l)^ZiZ2 

'  <  0 

16rjr2 


(52) 


as  lf(  0.  Therefore,  any  growing  mode  will  neutralize 
as  the  deformation  radius  drops  below  a  finite  thresh¬ 
old.  Such  stabilization  occurs  in  part  because  the  VR 
waves  on  opposite  edges  of  the  ring  become  shielded 
from  each  other,  and  can  no  longer  interact  effectively. 
Of  course,  for  more  subtle  reasons  the  growth  rate  may 
vary  nonmonotonically  with  decreasing  If,  before  shoot¬ 
ing  to  zero. 

Figure  3  illustrates  the  typical  properties  of  an  un¬ 
stable  mode  of  a  shallow-water  PV  ring.  The  angular 
phase  velocity  of  the  mode  is  greater  than  fij  but  less 
than  0,2-  Therefore,  the  inner  VR  wave  is  prograde, 
whereas  the  outer  VR  wave  is  retrograde.  As  predicted. 


^  To  begin  with,  all  modified  Bessel  functions  in  the  definition 
of  c  were  replaced  with  the  first  seven  terms  of  their  large- 
argument  (small  Ijj)  asymptotic  expansions  (Abramowitz  and  Ste- 
gun  1972).  The  resulting  approximation  for  c  was  then  reduced  to 
a  manageable  size  with  the  “Simplify”  routine  of  Mathematica  5.0. 
Equation  (52)  is  the  leading-order  term  of  the  approximation. 


0.2 
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0 
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Fig.  3.  The  dependence  of  the  complex  frequency  of  the  n  = 
3  growing  mode  of  a  QG  shallow-water  ring,  with  =  3/4,  r2  = 
5/4,  Z2  =  Zj  =  1,  and  thus  Zg  =  0.  Because  the  phase  velocity  Vfln 
is  between  Cl^  and  flj-  the  inner  VR  wave  is  prograde,  whereas  the 
outer  VR  wave  is  retrograde.  Note  also  that  the  growth  rate  Vj  is 
zero  for  all  Ijf  below  a  finite  threshold. 


the  mode  neutralizes  as  the  gravitational  restoring  force 
(parameterized  by  Ij^)  drops  below  a  finite  threshold.^ 

b.  A  category  3  hurricane 

We  now  reexamine  the  eyewall  instability  of  the  cat¬ 
egory  3  “hurricane”  in  Nolan  and  Montgomery  (2002, 
hereafter  NM02).  The  reader  may  consult  NM02  for  an 
analytical  description  of  the  basic-state  cyclonic  vortic- 
ity  ^(r,  z).  The  vorticity  is  easily  integrated  for  azi¬ 
muthal  velocity  v{r,  z).  The  equilibrium  profiles  of 
pressure  and  density,  p(r,  z)  and  p(r,  z),  are  found  by  an 
iterative  algorithm  that  solves  hydrostatic  and  gradient 
wind  balance  simultaneously.  For  this  study,  a  Jordan 
mean  sounding  gives  the  thermal  boundary  condition  at 
r  =  200  km  (Jordan  1958).  In  addition,  the  Coriolis 
parameter  is  set  to  /=  5  X  10“^  s“^.  Figure  4  portrays 
the  basic  state  of  the  vortex  on  the  r-p  plane.  Note  that 
the  eyewall  (where  v  peaks)  has  a  ringlike  PV  distribu¬ 
tion  of  finite  depth.'^ 


^  By  very  similar  analysis,  it  is  readily  shown  that,  in  the  context 
of  QG  theory,  the  internal  VR  modes  of  a  density-stratified  baro- 
tropic  PV  ring  ultimately  stabilize  as  the  buoyancy  frequency 
tends  toward  zero. 

'*  “Cloud-resolving”  numerical  simulations  suggest  that  the 
midlevel  PV  of  a  mature  hurricane  is  more  broadly  distributed  in 
radius  (e.g.,  Chen  and  Yau  2001).  So,  in  the  middle  troposphere, 
the  sharply  peaked  PV  profile  of  the  NM02  hurricane  might  be 
atypical.  Observations  of  Isabel  2003  further  show  that  the  PV 
distribution  of  a  real  hurricane  can  have  multiple  maxima  (Bell 
and  Montgomery  2006,  manuscript  submitted  to  Mon.  Wea.  Rev.). 
We  leave  this  more  complex  situation  for  future  consideration. 
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Fig.  4.  Basic  state  of  the  Nolan-Montgomery  category  3  hurricane.  Contours  ranging  from  380.23  to  466.46  K 
were  removed  from  the  stratosphere  of  the  0p  plot  due  to  excessive  bunching.  In  all  plots,  the  radius  and  pressure 
axes  are  scaled  linearly. 


According  to  NM02,  if  the  vortex  is  dry,  its  eyewall  is 
unstable.  Moreover,  the  azimuthal  wavenumber  of  the 
fastest  growing  eigenmode  is  «  =  3.  In  order  to  under¬ 
stand  how  moisture  influences  the  fastest  growing 
mode,  we  here  consider  a  set  of  numerical  experiments 
with  variable  Y.  Specifically,  we  let 

Y(r,p)  =  1  -  8[1  -  Y„,(cp)]-  (53) 

As  the  single  parameter  b  increases  from  zero  to  one, 
the  vortex  goes  from  a  dry  state  to  a  moist  state  that  is 
neutral  to  slantwise  convection.  The  lower  right  panel 
of  Fig.  4  shows  Y„,  for  the  vortex  under  consideration. 


In  each  experiment,  the  vortex  is  initialized  with  an 
arbitrary  n  =  3  perturbation.  Appendix  C  describes  the 
numerical  model  that  is  used  to  evolve  the  disturbance 
forward  in  time.  In  all  cases,  a  growing  eigenmode  pre¬ 
vails  after  a  brief  adjustment  period.  To  be  clear,  the 
eigenmode  is  a  disturbance  of  the  form 

<j)' = -b  C.C.,  (54) 

and  likewise  for  all  other  perturbation  fields.  As  before, 
V  =  vj(  +  ivj,  in  which  Vf,  and  vj  are  real.  Figure  5 
illustrates  the  exponential  growth  of  the  mode  for  in- 
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time  (lO'^s) 

Fig.  5.  Domain  maximum  of  the  n  =  3  geopotential  perturba¬ 
tion  of  the  NM02  hurricane,  with  variable  moisture.  After  a  brief 
transition  period,  the  perturbation  grows  like  e"''.  The  values  of  Vj 
are  given  to  the  right  of  each  curve.  As  additional  moisture  (in¬ 
creasing  e)  brings  the  vortex  toward  slantwise  convective  neutral¬ 
ity,  the  instability  slows  down.  In  this  plot,  !<()'  is  normalized 
to  its  initial  value. 

creasing  levels  of  moisture.  Evidently,  cloud  coverage 
tends  to  neutralize  the  mode.^ 

We  note  that  the  growth  rate  of  the  dry  mode  (3.3  X 
10“"^  s“^)  exceeds  the  value  of  2.07  X  10“^^  s“^  that 
appears  in  NM02.  During  the  review  process,  D.  Nolan 
(2006,  personal  communication)  was  able  to  repeat  the 
computation  of  NM02  at  a  finer  spatial  resolution, 
which  is  now  practical.  With  more  than  double  the 
number  of  radial  grid  points,  the  growth  rate  of  the  dry 
mode  increased  to  2.7  X  10“"^  s“^.  The  remaining  dis¬ 
crepancy  between  this  value  and  our  result  is  likely 
caused  by  (i)  our  different  surface  boundary  condition 
(w'  =  0  on  a  constant  p  surface  as  opposed  to  no  stress 
on  a  constant  z  surface),  (ii)  our  use  of  the  hydrostatic 
approximation,  and  (iii)  the  absence  of  explicit  eddy 
diffusion  in  our  model. 

c.  Comparison 

In  many  ways,  the  eyewall  instability  of  the  NM02 
hurricane  resembles  the  instability  of  a  QG  shallow- 
water  ring.  To  begin  with,  whether  dry  or  moist,  the 
unstable  eyewall  mode  of  the  NM02  hurricane  appears 


to  comprise  counter  propagating  VR  waves.®  To  see 
this,  let  us  examine  its  structure.  The  time-asymptotic 
PV  perturbation  approximately  has  the  form 


n' =  n(r,/7)e‘'<”‘"“’"> -b  c.c.  (55) 

The  two  plots  in  Fig.  6  show  |n|  for  e  =  0  and  for  e  = 
0.8875.  For  8  =  0,  the  dry  case,  the  radial  variation  of 
|n|  has  two  distinct  peaks.  These  peaks  occur  at  the 
maximum  positive  and  negative  values  of  Sll/Sr.  For  e 
=  0.8875, 1  n|  has  a  similar  structure,  but  the  inner  wave 
is  fractured. 

In  both  plots,  the  dashed  arc  is  where  the  Doppler- 
shifted  mode  frequency. 


o-R(r,  P)  =  Vr  -  nCl,  (56) 

is  zero.  Above  the  arc  a-f,  >  0,  and  below  the  arc  cr^  < 
0.  Whether  dry  or  moist,  the  inner  wave  is  above  the 
arc,  and  is  therefore  prograde  relative  to  the  mean  flow. 
In  contrast,  the  outer  wave  is  underneath  the  arc,  and  is 
therefore  retrograde.  These  properties  are  consistent 
with  the  azimuthal  propagation  rule  for  VR  waves  in  a 
cyclone:  the  propagation  is  prograde/retrograde  if  Sll/Sr 
is  positive/negative  (Montgomery  and  Kallenbach 
1997). 

True  VR  waves  should  also  be  slow,  in  some  sense, 
relative  to  IG  oscillations.  One  quantitative  condition 
for  slowness  is  given  below: 


{^dv/dp)^  p  / pQ\^d/cpd  _ 

Td^/dpRAp) 


(57) 


The  above  condition  states  that  the  Doppler-shifted 
mode  frequency  should  be  much  less  than  the  minimum 
frequency  of  slantwise  IG  oscillations  (the  square  root 
of  the  denominator;  see  appendix  B).  The  contours  in 
Fig.  6  show  the  r-p  variation  of  for  the  dry  and  moist 
modes.  In  the  eyewall,  both  modes  appear  to  margin¬ 
ally  satisfy  condition  (57).^ 

We  note  that  if  ©^  <C  1,  terms  of  order  ©^  and  higher 
are  subdominant  in  the  perturbation  equations.  Ne¬ 
glecting  these  terms  leads  to  the  zero-order  asymmetric 
balance  model  (Shapiro  and  Montgomery  1993;  Mont¬ 
gomery  and  Shapiro  1995;  Montgomery  and  Franklin 
1998;  Ren  1999;  McWilliams  et  al.  2003).  Like  QG 
theory,  the  asymmetric  balance  model  filters  out  iner¬ 
tia-gravity  waves.  Unlike  QG  theory,  the  basic  state 
can  have  a  Rossby  number  of  order  unity  or  greater. 


^  Dunkerton  (1990)  examined  a  similar  divergent  barotropic 
instability  of  a  zonal  shear  flow.  By  comparison,  he  showed  that 
the  growth  rate  of  the  instability  can  first  increase  and  then  de¬ 
crease  as  the  static  stability  of  the  atmosphere  tends  toward  zero 
by,  say,  increasing  cloudiness. 


^  Nolan  and  Montgomery  (2002)  discuss  the  Rossby  wave  na¬ 
ture  of  the  instability  for  the  dry  case.  The  objective  here  is  to 
extend  the  discussion  to  cases  in  which  clouds  are  present. 

’  We  should  point  out  that  as  T  — >  becomes  infinite 

everywhere  but  on  the  contour  cTp  =  0.  Therefore,  the  area  in 
which  condition  (57)  holds  must  vanish  as  s  — >  1. 
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Fig.  6.  The  PV  amplitude  |n|  (normalized  to  its  maximum)  and  slowness  parameter  (contours)  for  the  most  unstable  n  =  3 
eyewall  modes  of  a  dry  and  a  moist  NM02  hurricane.  In  both  plots,  the  radius  and  pressure  axes  are  scaled  linearly. 


In  summary,  the  pertinent  instabilities  of  a  QG  shal¬ 
low-water  ring  and  of  the  NM02  hurricane  both  consist 
of  slow,  counter  propagating  VR  waves.  By  analogy  to 
the  QG  shallow-water  ring,  it  is  no  surprise  that  reduc¬ 
tion  of  the  gravitational/buoyancy  restoring  force  by 
adding  moisture  suppresses  the  NM02  hurricane  insta¬ 
bility. 

Of  course,  the  analogy  has  limits.  To  begin  with,  the 
NM02  hurricane  mode  is  not  globally  slow,  as  it  is  in 
QG  theory.  Figure  6  indicates  severe  violation  of  con¬ 
dition  (57)  in  the  periphery  of  the  storm,  in  the  eye  of 
the  storm,  and  in  the  upper  troposphere.  Previous  stud¬ 
ies  have  shown  that  such  external  violation  can,  in  some 
cases,  notably  affect  the  growth  rate  of  a  vortex  mode 
(Schecter  and  Montgomery  2003).  Furthermore,  mois¬ 
ture  significantly  alters  the  r-p  structure  of  the  mode. 
Figure  6  alone  illustrates  a  vertical  elevation  and  scale 


reduction  of  PV  peaks.  Figure  7  illustrates  the  changes 
of  other  perturbation  fields.  Note  in  particular  that 
moisture  causes  mesoscale  spatial  oscillations  to 
emerge  in  the  eye.  The  emergence  is  most  striking  for 
a,'.» 

7.  The  eigenmode  problem  for  a  cloudy  vortex 

In  section  8,  we  will  examine  the  discrete  eigenmodes 
of  a  cloudy  vortex  that  resembles  an  incipient  tropical 
cyclone.  Before  doing  so,  it  is  helpful  to  discuss  the 
eigenmode  problem  more  generally. 

a.  Eigenmode  problem  for  a  generic  vortex 

Substituting  a  traveling  wave  ansatz  [Eq.  (54)  and 
relatives]  into  the  linearized  equations  of  motion  [Eqs. 
(29)-(34)],  and  eliminating  all  perturbation  fields  but 
the  geopotential,  we  obtain 


r  dr  \T  dp  dr  / 


d 

^  dp 


2d  dp  d^<P 


dB„ 


1 


r^T 


+  _ 

dp 


dp  \  r  dpdr 


1  do 

111 

(  r  dOp) 

,  dd* 

dp  ^ 

[r  dr) 

'  dr 

r  dr 

[t  dr) 

rtr  L  \  r  dr  )  \  F  dp  )  _ 


$  =  0, 


(58) 


in  which 


cr  =  V  —  nil, 


and 


(59) 


-  S0p 
dp 


/  5u\^  p  / pQ\^d/‘:pd 

Vdpj  Ra\p) 


(61) 


is  the  complex  Doppler-shifted  mode  frequency. 


®  A  more  realistic  representation  of  a  hurricane  would  have  a 
relatively  dry  eye  compared  to  the  eyewall,  in  which  case  the 
striking  fluctuations  of  oi'  at  small  radii  would  likely  disappear. 


(60) 
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e  =  0.0  0.75  0.8875 


density  potential  temperature  <  6'  ) 


azimuthal  velocity  (  v' ) 


Fig.  7.  Vertical  slices  of  the  unstable  n  =  3  eyewall  mode  of  the  NM02  hurricane,  for  increasing  levels 
of  moisture  from  left  to  right.  The  radius  and  pressure  axes  are  scaled  linearly,  and  have  the  same  limits 
for  every  plot  (see  top  left).  The  azimuth  (cp)  for  each  column  is  where  0^  has  its  peak  value.  The  wave 
amplitudes  are  arbitrarily  normalized  so  that  the  maximum  of  |  (()'  | ,  at  the  azimuth  shown,  is  100  m^  s“^. 
Bright  (dark)  shades  represent  positive  (negative)  values  of  the  plotted  fields.  The  dashed  curves  are  the 
zero  contours.  Because  of  excessive  bunching,  negative  contours  were  removed  from  the  e  =  0.8875 
velocity  plots. 
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The  values  of  v  that  permit  nontrivial  solutions  of  Eq. 
(58),  that  is,  the  eigenfrequencies,  are  constrained  by 
the  radial  and  vertical  (pressure)  boundary  conditions. 
We  will  specify  these  conditions  when  the  need  arises. 

b.  Eigenmode  problem  for  a  barolropic  vortex  with 
separable  T 


Clearly,  decreases  as  the  vertical  wavenumber  or  the 
Coriolis  parameter  increases;  also  varies  radially. 

Suppose  that  S’  =  cdOpIdp,  in  which  c  is  a  constant. 
Then, 


Equation  (58)  is  complicated,  and  useful  only  if  sim¬ 
plified.  For  example,  suppose  that  the  basic  state  of  the 
vortex  is  barotropic;  that  is. 


dv  30 
dp  dr 


(62) 


Then,  the  eigenmode  equation  reduces  to  the  following: 
1  d  (  r  SchX  1  d  / 


r  dr  \  ip  dr  J  ra  dr  \  ip 


r^ip 


(63) 


in  which 


P  / 


dp 


(64) 


RfCdQp/dp 

Further  suppose  that  the  buoyancy  reduction  factor 
is  separable.  Then,  we  may  write 


T^=%ir)9(p).  (65) 

Equation  (65)  permits  us  to  seek  eigenfunctions  of  the 
form 


<1>  =  (66) 

Upon  substituting  Eq.  (66)  into  Eq.  (63),  we  find  that  a 
nontrivial  solution  for  cp;^  can  exist  only  if 


r  P  ( PoY‘''‘'‘’'‘ K 

dp  iRd'Pip)  \  P  /  dp 


(67) 


in  which  A  is  a  constant.  The  vertical  boundary  condi¬ 
tions  determine  the  permissible  values  of  A.  Fligher  val¬ 
ues  of  A  correspond  to  eigenfunctions  with  smaller  ver¬ 
tical  wavelengths.  If,  for  example,  w'  =  0  at  p  = 
and  at  p  =  then  the  boundary  conditions  are  d'^J 
dp  =  Oat  p  =  p^i^  and  at  p  =  p^^^. 

Inserting  Eqs.  (65)-(67)  into  Eq.  (63),  and  dividing 
through  by  yields 


I  d  /  r 
r  dr  \  ijj  dr 


rcrdryip 


-  - 

^  Hf 


$,  =  0. 


(68) 


The  symbol  Ij^  here  denotes  the  internal  Rossby  defor¬ 
mation  radius,  defined  by 


Ir  =  Ym  IP,  (70) 

in  which  is  the  internal  Rossby  deformation  radius 
in  the  absence  of  moisture.  Since  Y  <  1,  it  is  clear  that 
moisture  reduces 

To  complete  the  eigenmode  problem,  we  must 
specify  two  radial  boundary  conditions.  At  the  origin, 
we  require  that  is  zero  (n  >  1)  or  finite  (n  =  0).  As 
r  ^  oo,  we  impose  either  a  radiation  or  evanescence 
condition,  depending  on  whether  the  mode  is  fast  or 
slow  relative  to  an  ambient  inertial  oscillation.  In  ap¬ 
pendix  D,  we  give  both  the  inner  and  outer  asymptotic 
forms  of  From  these,  we  can  infer  relations  between 
cp;^  and  dt^Jdr  near  the  origin  and  at  distant  radii.  We 
use  these  relations  for  mixed  boundary  conditions  in 
our  numerical  eigenmode  solver,  which  is  based  on  the 
shooting  method.  The  computational  boundaries  are 
typically  at  O.Olro  and  15^0,  in  which  r^  is  the  radius  of 
the  vortex  core. 

Interested  readers  may  consult  appendix  E  for  a  set 
of  polarization  equations  that  give  the  wavefunctions 
for  m',  v' ,  and  dp  in  terms  of  <l)x(0  ’F;^(p). 

8.  Suppression  of  spiral  inertia-gravity  wave 
radiation  and  improvement  of  vortex  resilience 
by  cloud  coverage 

During  their  early  stages  of  development,  and  per¬ 
haps  after  eyewall  breakdown,  tropical  cyclones  are 
presently  thought  to  have  PV  distributions  that, 
roughly  speaking,  decrease  monotonically  away  from 
the  storm  center.  In  this  section,  we  examine  how  mois¬ 
ture  affects  the  propagation  and  exponential  growth  of 
discrete  VR  waves  in  an  ageostrophic  monotonic  cy¬ 
clone.  The  most  familiar  of  these  waves  include  rotating 
tilts  and  elliptical  deformations  (e.g.,  Kelvin  1880; 
Schecter  and  Montgomery  2004).  They  are  distin¬ 
guished  from  continuum  VR  waves  (e.g.,  Montgomery 
and  Kallenbach  1997;  Bassom  and  Gilbert  1998)  in  that 
their  PV  fields  do  not  characteristically  suffer  spiral 
windup.  Typically,  spiral  windup  of  PV  leads  to  alge¬ 
braic  decay  of  the  asymmetric  geopotential  perturba¬ 
tion,  as  opposed  to  exponential  decay  or  growth. 

Discrete  VR  waves  have  some  role  in  regulating 
tropical  cyclone  intensity.  This  role  is  not  completely 
understood,  but  some  basic  theory  has  been  developed. 
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Fig.  8.  A  tanh  barotropic  cyclone  with  variable  cloudiness,  (a)  Relative  vorticity  f  =  tj  —  /, 
and  buoyancy  reduction  factor  T  for  increasing  levels  of  moisture  (increasing  e)  near  the 
radius  of  maximum  wind,  (b)  The  outer  relative  vorticity  has  a  small  but  finite  negative 
gradient.  In  both  (a)  and  (b),  f  is  normalized  to  its  central  value,  and  r  is  normalized  to  the 
core  radius  r„. 


In  linear  theory,  discrete  VR  waves  are  compelled  to 
decay  by  resonant  transfer  of  wave  activity  into  a  criti¬ 
cal  layer  (Briggs  et  al.  1970;  Pillai  and  Gould  1994; 
Schecter  et  al.  2000,  2002;  Balmforth  et  al.  2001; 
Schecter  and  Montgomery  2004,  2006;  Mallen  et  al. 
2005).  During  decay,  the  wave-mean  flow  interaction 
usually  increases  the  maximum  tangential  wind  speed 
in  the  vortex  core  (e.g.,  Montgomery  and  Enagonio 
1998;  Reasor  and  Montgomery  2001).  On  the  other 
hand,  if  >  /,  discrete  VR  waves  are  compelled  to 
grow  by  exciting  outward  propagating  spiral  IG  waves 
in  the  environment  (Ford  1994a, b;  Plougonven  and 
Zeitlin  2002;  Schecter  and  Montgomery  2004,  2006). 
This  spontaneous  radiation  transports  energy  and  an¬ 
gular  momentum  away  from  the  vortex.  The  following 
demonstrates  that  core  moisture  increases  the  absorp¬ 
tion  rate  of  wave  activity  in  the  critical  layer  (indirectly) 
and  thereby  inhibits  spontaneous  IG  wave  radiation.  In 
this  sense,  moisture  can  fortify  the  cyclone’s  intensity. 

The  following  demonstration  also  implies  that  cloud 
coverage  improves  the  resilience  of  a  monotonic  cy¬ 
clone  that  is  exposed  to  ambient  vertical  shear.  Suppose 
that  the  discrete  VR  wave  under  consideration  is  a  tilt 
mode  (n  =  1).  As  mentioned  previously,  adding  clouds 
to  the  cyclone  either  causes  or  accelerates  decay  of  that 
mode.  Therefore,  a  cloudy  monotonic  cyclone  realigns 
more  rapidly  than  its  dry  counterpart  after  it  is  tilted  by 
transient  vertical  shear.  As  we  will  see,  cloud  coverage 
also  boosts  the  azimuthal  phase  velocity  of  the  mode. 
Intuition  and  basic  modeling  suggest  that  sustained  ver¬ 


tical  shear  becomes  relatively  harmless  if  the  phase  ve¬ 
locity  of  the  mode  becomes  greater  than  the  shear  rate 
(e.g.,  Jones  1995;  Reasor  et  al.  2004). 

a.  A  barotropic  cyclone 

We  focus  the  present  study  on  the  discrete  VR  waves 
of  a  barotropic  cyclone,  with 


ar)  =  2/ 


1  —  tanh 


(71) 


As  r  increases  from  the  origin,  |  stays  roughly  constant, 
until  dropping  rapidly  toward  zero  near  r^  (see  Fig.  8). 
The  width  of  the  transition  layer  is  of  order  d  = 
0.0375^0.  The  Rossby  number  of  the  cyclone,  Ro  = 
Ci0)/f  =  4,  is  characteristic  of  an  incipient  storm.  Most 
importantly,  it  exceeds  unity.  The  angular  rotation  fre¬ 
quency  fi  (not  shown)  is  roughly  constant  for  r  <  r^, 
and  decays  nearly  as  r~^  for  r  >  r„. 

We  further  assume  that  the  buoyancy  reduction  fac¬ 
tor  does  not  vary  with  p.  Specifically,  we  let 

Y(/-)  =  [1  -  (72) 

which  has  the  minimum  value  (1  —  ef'  at  r^,  and  is  unity 
for  \  r  —  rg\:^2d  (see  Fig.  8).  The  8  parameter  can  have 
any  value  between  zero  and  one.  Note  that  our  choice 
of  Y  is  minimal  in  the  region  of  maximum  azimuthal 
wind  speed.  This  property  was  motivated  by  the  simu¬ 
lated  hurricane  data  in  Fig.  2. 


328 


JOURNAL  OF  THE  ATMOSPHERIC  SCIENCES 


Volume  64 


b.  The  growth  rate  of  a  discrete  vortex  Rossby  wave 

Figure  9  shows  a  horizontal  slice  of  the  geopotential 
perturbation  f)'  of  a  discrete  VR  wave  in  a  cloudy  tanh 
cyclone  [Eqs.  (71)-(72)].  The  oscillation  frequency 
of  the  wave  exceeds  the  Coriolis  parameter  /.  Conse¬ 
quently,  the  VR  wave  resonantly  excites  an  outward- 
propagating  spiral  IG  wave  in  the  environment.  The 
coupled  VR  and  IG  waves  constitute  a  single  eigen¬ 
mode. 

It  has  been  shown  that,  near  marginal  stability,  the 
growth  rate  of  a  discrete  VR  wave  in  a  dry  barotropic 
vortex  is  given  by  the  formula  (cf.  Schecter  and  Mont¬ 
gomery  2004) 

V/  =  I'rad  -  (73) 

in  which 

rl 

I^rad-^9^[t7xn],  (74) 

and 

TT  r^\Uj^\^dr\/dr 

M  _  \nd£l/dr\ 

ITere,  UJj)  and  V^{r)  are  the  radial  wavefunctions  of 
the  radial  and  azimuthal  velocities  [see  appendix  E  for 
their  connection  to  ‘i5;^(r)].  The  superscript  *  is  the  com¬ 
plex  conjugate  operator,  and  9i[. . .]  is  the  real  part  of 
the  quantity  in  square  brackets.  The  critical  radius  r^  is 
defined  by 

^^(g)  -  (76) 


It  is  where  the  angular  rotation  frequency  of  the  vortex 
equals  the  angular  phase  velocity  of  the  wave.  The  ra¬ 
dius  r„  defines  the  outer  vortex  boundary.  Its  exact 
value  is  somewhat  arbitrary,  but  must  exceed  r^.  The 
denominator  M  is  defined  by 


M=  - 


dr  |n-|2 


—  -H  C.C.) 


(77) 


in  which  -f  denotes  a  radial  integral  that  excludes  an 
“infinitesimal”  band  at  r^.  Provided  that  our  choice  of 
is  not  orders  of  magnitude  greater  than  r^,  M  is  gen¬ 
erally  positive  (Schecter  and  Montgomery  2004,  2006). 
Vertical  boundaries  do  not  appear  explicitly  in  the 
equations  above.  However,  the  hidden  vertical  bound¬ 
ary  conditions  are  w'  =  0  (and  hence  0p  =  0)  at  p  = 

Pmax  Hud  at  p  Pmin- 

Equation  (73)  indicates  that  the  growth  rate  of  a  dis¬ 
crete  VR  wave  has  two  parts.  The  first  term,  is 


Fig.  9.  Horizontal  slice  of  the  geopotential  perturbation  4>'  of 
the  n  =  2  discrete  VR  wave  of  a  tanh  barotropic  cyclone  (Fig.  8), 
and  the  IG  wave  radiation  that  it  generates.  The  moisture  param¬ 
eter  for  this  case  is  s  =  0.725,  the  Rossby  number  is  Ro  =  4,  and 
the  dry  Rossby  deformation  radius  is  l'^  =  r„.  Solid/dashed  con¬ 
tours  correspond  to  positive/negative  values  of  (j)'-  The  contour 
levels  are  ±0.05,  ±0.3,  ±0.75,  and  ±1.5,  in  units  that  are  normal¬ 
ized  to  one-half  the  amplitude  of  cj>'  at  r  =  /■„. 


proportional  to  the  outward  angular  momentum  flux  of 
the  spiral  IG  wave  that  the  discrete  VR  wave  emits  into 
the  environment.  This  flux  is  generally  positive;  there¬ 
fore,  the  IG  wave  emission  destabilizes  the  discrete  VR 
wave.  The  second  term,  —  v^i,  accounts  for  the  resonant 
feedback  of  PV  stirring  in  the  critical  layer.  It  is  nega¬ 
tive,  since  drjidr  =  dj/dr  <  0,  and  therefore  acts  to 
stabilize  the  discrete  VR  wave. 

The  preceding  growth  rate  formula  also  applies  to 
the  discrete  VR  waves  of  a  cloudy  barotropic  vortex  in 
which  Y  is  separable  [see  Eq.  (65)].  Moisture  implicitly 
changes  the  right-hand  side  of  Eq.  (73)  by  altering  the 
wavefunctions,  moving  r^,  and  giving  Ifi  radial  variation 
[see  Eq.  (70)].  Readers  who  would  like  to  verify  Eq. 
(73)  in  the  moist  context  may  follow  the  dry  derivation 
in  Schecter  and  Montgomery  (2004).  In  doing  so,  one 
must  replace  the  dry  Boussinesq  wave  activity  equation 
with  its  moist  non-Boussinesq  counterpart.  Appendix  F 
provides  the  pertinent  equation. 

c.  Computational  results 

Figures  10  and  11  illustrate  the  effects  of  moisture  on 
the  discrete  VR  waves  of  the  tanh  cyclone  that  is  under 
consideration  [Eqs.  (71)  and  (72)].  These  waves  corre¬ 
spond  to  (nonsingular)  numerical  solutions  of  the 
eigenmode  problem  of  section  7b.  Figure  10  shows  how 
the  critical  radii  vary  with  the  moisture  parameter  e. 
The  four  plots  cover  azimuthal  wavenumbers  IM.  The 
different  curves  in  each  plot  correspond  to  different 
values  of  the  dry  Rossby  deformation  radius,  In 
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Fig.  10.  Critical  radius  vs  the  moisture  parameter  e  for  a  spectrum  of  discrete  VR  waves  of  a  tanh 
barotropic  cyclone,  with  Ro  =  4.  In  general,  decreases  toward  the  core  radius  /•„  of  the  cyclone  as  the 
level  of  moisture  increases.  All  lengths  are  normalized  to  r„. 


general,  adding  clouds  increases  the  azimuthal  phase 
velocity  of  a  wave  and  thereby  decreases  toward  i-q.® 

Figure  11  shows  how  the  growth  rates  vary  with  8.  In 
general,  adding  clouds  tends  to  decrease  vj  and  stabilize 
a  wave.  Stabilization  occurs  because  cloud  coverage 
brings  the  critical  radius  to  a  region  of  relatively  large 
negative  PV  gradient  {(Fqldr).  This  enables  critical  layer 
damping  to  prevail  over  radiative  pumping  (see  below). 

Figure  12  verifies  the  analytical  growth  rate  formula 
of  section  8b  for  a  particular  VR  wave  that  is  near 
marginal  stability.  Again,  the  growth  rate  is  plotted  as  a 
function  of  cloudiness  in  the  cyclone.  The  X  symbols 


®  Note  that  increasing  the  phase  velocity  of  a  discrete  VR  wave 
corresponds  to  slowing  its  retrograde  propagation  relative  to  the 
mean  flow  in  the  vortex  core. 


correspond  to  Eq.  (73).  To  evaluate  the  right-hand  side 
of  this  equation,  we  used  the  wavefunctions  {U^,  V;^, 
ch;^)  and  critical  radius  (r^)  of  the  computed  eigenmode. 
To  evaluate  M,  we  removed  the  interval  \r  -  rj  < 
2Q\vilndVtldr^_\  from  the  integral;  that  is,  we  removed 
an  estimate  of  the  linear  critical  layer  (cf.  Schecter  et  al. 
2000;  Schecter  and  Montgomery  2004).  Evidently,  Eq. 
(73)  accurately  reproduces  the  values  of  i'/  (+  symbols) 
that  are  obtained  directly  from  the  eigenmode  solu¬ 
tions. 

Figure  12  also  shows  how  the  radiative  pumping  rate 
(Aad)  critical  layer  damping  rate  (— Vcz)  vary  with  the 
moisture  parameter  e.  To  obtain  these  rates,  and  v^i 
were  evaluated  with  As  e  increases  from  zero  to 

one,  the  radiative  pumping  rate  remains  roughly  con¬ 
stant.  However,  the  negative  PV  gradient  at  explosively 
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Fig.  11.  Growth  rate  Vj  [normalized  to  {(0)]  vs  the  moisture  parameter  e  for  a  spectrum  of  discrete  VR 
waves  of  a  tanh  barotropic  cyclone,  with  Ro  =  4.  In  all  cases,  Vj  eventually  becomes  negative  as  the  level 
of  moisture  increases.  All  deformation  radii  are  normalized  to  the  core  radius  r„. 


grows,  and  the  critical  layer  damping  rate  increases  by 
orders  of  magnitude.  Consequently,  vj  becomes  negative. 

We  note  that  as  the  growth  rate  becomes  negative, 
the  discrete  VR  wave  transforms  from  an  eigenmode 
into  a  quasimode.  Quasimodes  are  not  exact  solutions 
to  the  eigenmode  problem.  For  instance,  the  PV  distur¬ 
bance  of  a  quasimode  grows  with  time  in  the  thin  criti¬ 
cal  layer  centered  at  r^',  in  contrast,  the  PV  disturbance 
of  a  damped  eigenmode  would  decay  everywhere.  The 
numerical  method  for  finding  quasimodes  is  similar  to 
the  numerical  method  for  finding  eigenmodes,  and  has 
been  described  elsewhere  (e.g.,  Spencer  and  Rasband 
1997;  Schecter  and  Montgomery  2004,  2006). 

We  must  also  emphasize  that  the  above  discussion  is 
based  on  linear  theory.  Nonlinear  processes  can  revive 
a  linearly  damped  VR  wave.  For  instance,  dry  studies 
have  shown  that  radiative  pumping  ultimately  over¬ 
takes  critical  layer  damping  if  the  initial  VR  wave  am¬ 
plitude  exceeds  a  finite  threshold  (Schecter  and  Mont¬ 
gomery  2006). 


9.  Concluding  remarks 

In  this  paper,  we  presented  a  basic  model  for  small- 
amplitude  wave  dynamics  in  a  cloudy  vortex.  The  wave 
equations  were  derived  with  the  following  set  of  ap¬ 
proximations: 

(i)  an  air  parcel  conserves  its  total  water  content  (liq¬ 
uid  plus  vapor); 

(ii)  the  mixing  ratios  of  water  vapor  and  liquid  cloud 
in  an  air  parcel  instantaneously  adjust  to  mm[q„ 
q^.^]  and  max[0,  q,  -  q^^],  respectively; 

(hi)  an  air  parcel  conserves  its  total  entropy; 

(iv)  the  wave  motions  are  in  hydrostatic  balance; 

(v)  the  wave  motions  follow  linearized  equations; 

(vi)  the  basic-state  vortex  has  no  secondary  circulation; 
Boundary  layer  drag  and  eddy  viscosity  were  also  ig¬ 
nored.  Clearly,  this  model  was  not  designed  to  precisely 
simulate  atmospheric  vortices.  Rather,  it  was  designed 
as  a  tool  to  easily  explore  the  most  basic  effects  of  cloud 
coverage  on  vortex  wave  dynamics. 
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Fig.  12.  Verification  of  the  analytical  growth  rate  formula  for  a 
discrete  VR  wave  in  a  barotropic  tanh  cyclone,  with  Ro  =  4.  The 
-I-  symbols  correspond  to  numerical  solutions  of  the  eigenmode/ 
quasimode  problem,  and  the  X  symbols  correspond  to  Eq.  (73). 
The  dashed  curves  show  the  radiative  pumping  rate  and 

critical  layer  damping  rate  (  — Vci)-  Increasing  cloudiness  near  the 
radius  of  maximum  wind  (increasing  e)  ultimately  enables  critical 
layer  damping  to  prevail. 


The  moist  wave  equations  that  result  from  (i)-(vi) 
are  identical  to  the  linearized  dry  hydrostatic  primitive 
equations,  with  reduced  cooling/warming  by  updrafts/ 
downdrafts.  We  note  that  a  similar  model  was  recently 
reported  to  have  been  used  to  numerically  simulate  the 
effect  of  moisture  on  vortices  that  are  exposed  to  sus¬ 
tained  vertical  shear  (Patra  2004).  Here,  we  derived  the 
model  from  “first  principles,”  in  order  to  clarify  the 
underlying  assumptions.  In  addition,  we  provided  an 
explicit  formula  for  estimating  the  r-p  variation  of  the 
buoyancy  reduction  factor  Y  from  observations  [Eq. 
(32)]. 

The  latter  sections  of  this  paper  briefly  examined  the 
influence  of  moisture  on  eyewall  instability,  spontane¬ 
ous  IG  wave  radiation  and  vortex  resilience.  We  have 
by  no  means  closed  the  book  on  these  subjects.  Rather, 
we  have  presented  a  few  case  studies  to  motivate  future 
research.  In  one  case  study,  we  showed  that  cloud  cov¬ 
erage  can  temper  the  asymmetric  eyewall  instability  of 
a  hurricane-like  vortex.  In  another  case  study,  we 
showed  that  clouds  in  the  core  of  a  monotonic  cyclone 
can  cause  or  accelerate  the  decay  of  discrete  VR  waves. 
This  both  inhibits  the  spontaneous  radiation  of  IG 
waves  and  improves  vortex  resilience. 
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APPENDIX  A 


Thermodynamic  Formulas 

In  this  appendix,  we  derive  convenient  formulas  for 
the  moist  adiabatic  pressure  derivatives  of  density  po¬ 
tential  temperature. 


a.  Unsaturated  air 


For  unsaturated  air,  in  which  =  q„  we  may  write 


(Al) 


in  which  T“(p,  s^,  q,)  =  T{p,  q^  q„  q,).  It  follows 
that 


\  V  sp  A, 


(A2) 


We  must  now  evaluate  the  pressure  derivative  of  r“  at 
constant  s„^  and  q,. 

To  readily  exploit  a  published  result,  we  introduce 
the  variable  z(p,  s^,  q,),  defined  such  that 


^P  ' 

Using  Eq.  (A3),  we  may  write 


(A3) 


1 

/ar“ 

\  1  /ar“/ 

1 

'J-'U 

\  Sp  , 

~  gpT'^  V  / 

Rd^  +  1 

(  ST“' 

1 

1 

+ 

K  5z  , 

(A4) 


The  above  z  derivative  of  7“  is  minus  the  adiabatic 
lapse  rate  of  an  unsaturated  hydrostatic  atmosphere, 
with  constant  q,.  Standard  texts  [e.g.,  Emanuel  1994, 
Eq.  (4.7.4)]  derive  the  following  formula  for  this  lapse 
rate: 


g  l  +  q, 

\  Cpd  1  +  qApJCpa  ■ 


(A5) 


One  may  now  execute  the  following  chain  of  substi¬ 
tutions: 


(i)  substitute  Eq.  (A5)  into  Eq.  (A4); 

(ii)  substitute  Eq.  (A4)'  into  Eq.  (A2); 

(hi)  substitute  Eq.  (Al)  into  Eq.  (A2)'. 

Here,  primes  indicate  modified  equations.  The  last  sub¬ 
stitution  (iii)  yields  formula  (15)  of  the  main  text. 


b.  Saturated  air 

For  saturated  air,  in  which  q.^  =  q^^.{p,  s^,  q,),  we 
may  write 
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-J 


1 


gt 


We  must  now  evaluate  the  pressure  derivatives  of 
(A6)  and  at  constant  and  q^. 

To  begin  with,  we  note  that 


in  which  T^(p,  s„,  q,)  =  T{p,  s^,  q^  q^^,  q,).  It  follows 


that 

^  1 

dr\ 

^  1 

1 

1  ^ 

r' 

\  dp  j 

gpr  ' 

\  dZ  /, 

\  I  P  Cpop  T\  dp  ] 

\  /  s^.qt  L  \  / 


D  1  + 

_ 

gP  1  + 


'dq^ 


e  +  <7^, 


dp 


(A7)  in  which  z  is  defined  by  Eq.  (A3),  and 


(A8) 


g  1  + 

A„,9,  ^pd  1  +  (ly,*(^pJCpd 


1  + 


1  + 


p.Ci 


^pd  +  qvifCpv  RJT’‘)^{Cp^  +  q^^c  ^ 


(A9) 


The  above  formula  for  the  adiabatic  lapse  rate  of  a 
saturated  hydrostatic  atmosphere  is  derived  in  Eman¬ 
uel  [1994,  Eq.  (4.7.3)]. 

Next,  we  observe  that 


we  further  obtain 

(dq^^\^  L^q^^il  +  q^Je) 
\  dr  Ip  ~  K{rf 


(A14) 


e 


* 


(AlO) 


in  which  is  the  saturation  vapor  pressure.  Since,  to  an 
excellent  approximation,  is  a  function  of  tempera¬ 
ture  r  alone,  q^^  is  a  function  of  only  p  and  r.  Ac¬ 
cordingly,  we  may  write 

/dqi,^\  /  dq^^\ 

V  ^p  \  dp  Jr^  V  ar  /p\  dp 

/ar\ 

(All) 

Upon  differentiating  Eq.  (AlO)  with  respect  to  p,  we 
obtain 


dqr,:,\ 
dp  /  U 


+  9'„=K/e) 

p 


(A12) 


Upon  differentiating  Eq.  (AlO)  with  respect  to  r,  and 
using  the  Clausius-Clapeyron  equation. 


dr  ~  Rjirf  ’ 


(A13) 


One  may  now  execute  the  following  chain  of  substi¬ 
tutions: 

(i)  substitute  Eqs.  (A12)  and  (A14)  into  Eq.  (All); 

(ii)  substitute  Eq.  (A9)  into  Eqs.  (AS)  and  (All)'; 
(hi)  substitute  Eqs.  (AS)'  and  (All)"  into  Eq.  (A7); 
(iv)  substitute  Eq.  (A6)  into  Eq.  (A7)'. 

Here,  single  or  double  primes  indicate  equations  that 
have  been  modified  once  or  twice,  respectively.  Simpli¬ 
fying  Eq.  (A7)"  yields  formula  (16)  of  the  main  text. 

APPENDIX  B 


Symmetric  Stability 

In  the  local  axisymmetric  model  [Eqs.  (35)  and  (36)], 
the  perturbation  equations  (section  4b)  have  constant 
coefficients.  Accordingly,  a  general  disturbance  is  a  su¬ 
perposition  of  Fourier  modes,  which  have  the  form 
p^^-i(vt-kpp-k,r) ^  in  which  a  is  a  complex  constant.  The 
dispersion  relation  is 


ddpRd  (  P  \^d'Cpd 


^2=  - 


dp  P  \Po 


dv  _  _ 

-  2k  —  I  +  T)  I, 

dp 


(Bl) 


in  which  k  =  k^lkp.  Derivation  of  Eq.  (Bl)  employs  the 
thermal  wind  relation  [Eq.  (2S)]. 
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The  vortex  is  locally  stable  only  if  >  0  for  all  k. 
Considering  the  limits  |  k|  ^  oo  and  |  k|  ^  0,  we  obtain 
the  familiar  conditions  of  gravitational  and  inertial  sta¬ 
bility; 


dO 

Y — ^<0  and  ^|>0. 
dp 


(B2) 


Assuming  gravitational  stability,  is  minimized  when 


^dv/dp  p  ! pQ\^^dfCpd 
TdO/dp  Rd\pJ 

Substituting  Eq.  (B3)  into  Eq.  (Bl)  yields 


(B3) 


V 


2 

min 


{^dv/dpf  p  (pf\Rd/cpd 

TdO^/dp  rAp) 


(B4) 


Requiring  that  >  0  gives  the  condition  of  slantwise 
convective  stability;  that  is,  Eq.  (37)  of  the  main  text. 


APPENDIX  C 


Numerical  Simulation 

In  this  appendix,  we  briefly  describe  the  numerical 
model  that  is  used  to  integrate  the  linearized  wave 
equations  (29)-(34)  forward  in  time. 

a.  Prognostic  equations 

To  begin  with,  let  us  expand  each  perturbation  field 
in  an  azimuthal  Fourier  series.  Specifically,  let 


in  which  jB^,  g,  Sr^,  p^,  and  8p^  are  adjustable  constants. 
The  simulation  domain  is  <  p  <  Pmax  0  < 
r  <  r^. 


b.  Diagnostic  equations 

At  each  time  step,  we  must  solve  two  diagnostic 
equations  in  order  to  evaluate  the  right-hand  sides  of 
Eqs.  (C2)-(C4).  The  first, 

Pmin.  0  +  3„(r,  p,  t),  (C6) 

is  derived  from  hydrostatic  balance.  It  involves  the  fol¬ 
lowing  weighted  pressure  integral  of  the  density  poten¬ 
tial  temperature  perturbation: 


3„  = 

The  second. 


P  Y^/cpd 


(Cl) 


oj„  =  (o„(r,  p„i„,  t)  +  y„(r,  p,  t),  (C8) 

is  derived  from  mass  continuity.  It  involves  another 
pressure  integral. 


1  d(ru„)  in 

dp 

Pmin 

- 1 

1 

+ 

1 

_ 1 

(C9) 


Clearly,  the  right-hand  sides  of  Eqs.  (C6)  and  (C8)  also 
depend  on  the  vertical  boundary  conditions. 

For  simplicity,  we  here  let 


(j)'  =  Mr,p,  0  +  S  l<l^nir,P,  t)e‘'"^  +  c.c.],  (Cl) 

n  =  l 

and  likewise  for  all  other  fields.  Our  numerical  model  is 
based  on  the  following  prognostic  equations  for  the 
Fourier  coefficients: 

and 


(C2) 

dv  in 

-pUn  - 

(C3) 

(^n(r,  Pmini  t)  =  0, 


(CIO) 


and 


w„(r,  p„ax-  0  =  0.  (Cll) 

Constraint  (CIO)  leaves  only  on  the  right-hand  side 
of  Eq.  (C8).  Since  constraint  (Cll)  must  hold  at  all 
times,  we  require  that  at  p^ax  i®  initially  zero,  and 
that 


C  Pmax  ^ 

1  d{ru„)  in 

J,  . 

r  dr  r 

(C12) 


f  d  _  \  dO  dd 

Note  that  we  have  added  a  sponge,  with  Rayleigh 

damping  coefficient  j3(r,  p),  to  absorb  inertia-gravity 
wave  radiation  at  the  top  and  outer  radius  of  the  simu¬ 
lation  domain.  The  damping  coefficient  takes  the  form 


1  -H  tanh 


1  +  tanh 


We  may  use  the  linearized  velocity  equations,  (C2)  and 
(C3),  to  evaluate  the  time  derivative  that  appears  in  the 
integrand.  Upon  doing  so,  we  obtain 

f  Pmax  C 

Pmini  0  =  “  dp  - - i  (C13) 

Pmin  /'max  /'min 

in  which 


1  a  a  n^ 

r  dr  ^  dr  ’ 


(C5) 


(C14) 
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i  d  _  inrt  n  _ 

^~~rYr  ~  iji)v„ 


1  d 
r  dr 


in  dv 


(C15) 


Equation  (C13)  is  a  second-order  ordinary  differential 
equation  in  r  for  the  function  4>„{r,  Pmiin  0-  To  solve  it, 
we  must  specify  boundary  conditions  at  r  =  0  and  at 
r  =  r^. 


c.  Radial  boundary  conditions 

At  the  domain  radius  r^,,  which  is  far  from  the  vortex 
center,  we  impose  an  artificial  wall,  that  is,  we  set 

M„(r6,  p,  t)  =  0.  (C16) 

The  radial  velocity  equation  (C2)  automatically  main¬ 
tains  zero  radial  velocity  at  r^,  at  all  times  if  we  impose 
the  condition 

(/-ft,  p,  t)  =  ^{r,,,  p)vj^r^,  p,  t)  (C17) 


The  parameters  of  the  “hurricane”  simulations  of 
section  6  are  listed  below: 

N,  =  400;  ^r  =  500  m; 

Np  =  397;  Ap  =  230.66  Pa; 

Tft  =  200  km; 

=  6115.12  Pa;p„3,  =  97  457.47  Pa; 
r^  =  175  km;  Sr^  =  5  km; 

=  15  000  Pa;  8p^  =  2500  Pa; 
j3,  =  0.03  s“^  At  =  2.5  s. 

Further  details  of  the  discretization  are  too  lengthy  for 
discussion.  Suffice  it  to  say  that  the  results  of  section  6 
appear  to  be  converged  with  respect  to  increasing  spa¬ 
tial  and  temporal  resolutions.  In  addition,  the  results 
have  no  appreciable  change  upon  bringing  j3j  to  zero. 

The  numerical  algorithm  has  passed  a  series  of  tests, 
in  which  simulated  eigenmodes  and  quasimodes  were 
shown  to  behave  properly.  Specifically,  their  complex 
frequencies  matched  those  obtained  from  the  shooting 
method. 


on  the  geopotential  perturbation.  Flydrostatic  balance 
further  implies  that 

dQ  p  { PfA^^d/Cpd  d 

^(/-ft,p,0=  -[|(/-ft,pK(rft,p,0]. 

(CIS) 

We  impose  no  extra  condition  on  the  tangential  veloc¬ 
ity  v„  at  rft. 

At  the  origin  {r  0)  we  enforce  the  usual  constraints 
on  the  Fourier  components  of  the  perturbation  fields. 
For  the  scalar  fields,  we  have 


0p,«  ^0  for  n  >  1 

(C19) 

and 

d4>o  30p,o  ^ 

dr  ’  dr 

(C20) 

For  the  velocity  field,  we  have 

u„ ,  u„  ^  0  for  n  ¥=  1 

(C21) 

and 

^  iv-^. 

(C22) 

d.  Discretization  parameters 

The  computational  grid  has  radial  grid  points  and 
Np  pressure  points.  The  radial  grid  spacing  is  Ar,  and 
the  pressure  spacing  is  Ap.  The  first  radial  grid  point  is 
at  r  =  Ar  and  the  last  is  at  r,,.  The  first  pressure  grid 
point  is  at  p  =  Pmm,  and  the  last  is  at  p^ax-  The  equa¬ 
tions  are  marched  forward  in  time  with  a  fourth-order 
Runge-Kutta  scheme,  in  steps  of  length  Af. 


APPENDIX  D 


Asymptotic  Wavefunctions 

In  this  appendix,  we  present  the  asymptotic  (small 
and  large  r)  solutions  of  the  geopotential  wavefunction 
of  an  eigenmode  of  a  barotropic  vortex  that  has  a  sepa¬ 
rable  Y  distribution.  By  assumption,  the  modal  oscilla¬ 
tion  frequency  (v^)  is  positive. 

Near  the  origin,  the  only  regular  solution  of  the  geo¬ 
potential  wavefunction  is 


r^O.  (Dl) 

Recall  that  n  is  the  azimuthal  wavenumber. 

As  r  becomes  infinite,  there  are  two  possibilities.  Let 


v/f  -  1 


(D2) 


If  91[k2]  >  0,  then  the  modal  oscillation  can  resonate 
with  an  environmental  IG  wave.  Consequently,  the 
mode  is  radiative,  and 


<1>A  (D3) 

Flere,  y  =  kr,  —  tt/2  <  arg(y)  <  Tr/2,  and  is  a  Han- 
kel  function  of  the  first  kind. 

If  91[A:2j  <  0,  then  the  mode  is  presumably  evanes¬ 
cent.  Consequently, 

‘Ta  Koiy),  r  ^  00.  (D4) 

Here,  y  =  \/—k^r,  —tt/2  <  arg(y)  <  ttI2,  and  Kq  is  a 
modified  Bessel  function.  Note  that  the  n  dependence 
of  the  radial  eigenfunction  is  subdominant  at  large  ra¬ 
dii,  and  appears  in  neither  (D3)  nor  (D4). 
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APPENDIX  E 


For  a  barotropic  vortex,  one  can  readily  verify  that 


Polarization  Equations  for  an  Eigenmode  of  a 
Cloudy  Barotropic  Vortex  with  Separable  T 


/d  _  a 
\dt  d(p 


dr\ 
dr  ’ 


(F3) 


Solving  the  eigenmode  problem  of  section  7b  yields  a 
geopotential  wavefunction  of  the  form  ch  = 

The  wavefunctions  for  the  radial  and  azimuthal  veloc¬ 
ities  are  given  by 

i  (  d<l>>  \ 

Uir,p)  =  -  =  [.—  --^,)^^  (El) 

and 


1  /  ncr  dOv 

E(gp)  =  T  (  TT'I’x  - 


(E2) 


The  density  potential  temperature  wavefunction  is 


0(r,p) 


rAp)  ^  dp  ■ 


(E3) 


In  addition  to  the  principal  polarization  equations 
(E1)-(E3),  there  are  several  auxiliary  relations.  The 
pressure  (vertical)  velocity  wavefunction  is  given  by 


(E4) 


and  the  PV  wavefunction  by 


n(G  p) 


1  d{rV)  in 

- 7 - U 

r  dr  r 


de^  30 
^-T]  — .  (E5) 
dp  '  dp 


Here,  we  have  placed  an  aberrant  circumflex  over  the 
wavefunction  to  prevent  confusion  with  our  symbol  for 
total  PV,  n(r,  (f>,  p,  t).  Finally,  the  wavefunction  of  the 
pseudo  PV  perturbation,  defined  in  appendix  F,  is 


„  i  dj] 


APPENDIX  F 


Wave  Activity  in  a  Cloudy  Barotropic  Vortex 


Our  construction  of  a  wave  activity  principle  requires 
that  we  first  introduce  a  pseudo  potential  vorticity  per¬ 
turbation. 


n;- 


d  (  dd 
n'F,;^F0;r,-ln(T^ 


_  -  a  /  e’  \  , 

dp\Td0^/dpl 

in  which 


(FI) 


V 


1  d(rv')  1  du' 

r  dr  r  dcp  ' 


(F2) 


in  the  context  of  linear  theory. 

A  convenient  local  measure  of  wave  activity  is  the 
moist  angular  pseudomomentum, 

r{n'„f  r  dv' 

r  = - - - -  n'  - 

P  2d^/dr  Tddp/dp  P  dp-  ^  ^ 

The  moist  angular  pseudomomentum  satisfies  the  flux 
conservative  equation, 

^TV.ir,  =  0,  (F5) 


in  which 


!Fp—  —  ru  V  r  +  -j  vLp  +  r 


{u'f  {v’f 


\Poj  2pTde  /dp}^ 


rr} 


Tddp/dp 


u'd'p  +  rv'o)'  |z. 


(F6) 


The  reader  may  verify  Eq.  (F5)  in  a  straightforward 
manner,  using  Eqs.  (30),  (31),  (33),  (62),  and  (F3). 

Integrating  Eq.  (F5)  over  a  cylinder  of  radius 
bounded  above  and  below  at  Pniax>  yields 


dA  fPmax 
dt 


Pmin  *'  "n" 


dp  d(p\r^u’v'\ 
r^T) 


dr  d(p 


u'd'  +  r^u'co' 
Yddp/  dp  P 


Pmax 

? 

-*Pmin 

(F7) 


in  which 


A  ^ 


Pmax  r  r 


Pmin  J  0  J  TT 


dp  dr  d(prL  (F8) 


is  the  total  wave  activity  in  the  cylinder.  The  total  wave 
activity  equation  (F7)  is  crucial  for  deriving  the  growth 
rate  formula  [Eq.  (73)]  of  a  discrete  VR  wave  in  a 
barotropic  vortex,  following  the  procedure  of  Schecter 
and  Montgomery  (2004,  2006). 

Suppose  that  0p  =  0  at  p^^^  and  Pniax-  Then,  by  the 
linearized  heat  equation  (31)  of  a  barotropic  vortex,  w' 
would  also  be  zero  on  the  vertical  boundaries.  As  such, 
the  second  (radial)  integral  on  the  right-hand  side  of 
Eq.  (F7)  would  vanish. 
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APPENDIX  G 


The  Meaning  of  Y 

In  this  appendix,  we  clarify  the  motivation  for  calling 
Y  the  “buoyancy  reduction  factor.”  This  nomenclature 
alludes  to  the  reduction  of  the  vertical  (z)  buoyancy 
restoring  force  that  a  parcel  feels  after  a  moist  adiabatic 
vertical  displacement,  or  pressure  change  (cf.  Smith 
et  al.  2005).  The  vertical  buoyancy  acceleration  is  de¬ 
fined  by 

B^-g-,  (Gl) 

P 


in  which  p'  is  the  local  perturbation  from  the  basic-state 
density 


P 


Cpd 


(G2) 


Consider  a  moist  fluid  parcel  that  moves  from  the 
pressure  level  p  —  Sp  to  p.  We  may  write  the  density 
change  of  the  parcel  as  follows: 


Sp 


8p  + 


9p 

de. 


CpdP 


(G3) 


in  which  the  new  symbol  is  the  specific  heat  of  dry 
air  at  constant  volume,  and  the  coefficients  of  the  in¬ 
finitesimals  are  evaluated  at  the  original  parcel  posi¬ 
tion. 

Assume  that  the  parcel  initially  has  basic-state  values 
of  density  and  density  potential  temperature.  Further 
assume  that  the  parcel  keeps  its  radial  coordinate  r 
fixed.  Then,  the  density  perturbation  at  the  new  parcel 
position  is  given  by 


p'  =  p(r,p  -  Sp)  -h  Sp  -  p{r,p) 

c^P 


CpdP  Sp 


dp\  p 

Sp  --SPp. 


(G4) 


It  is  straightforward  to  show  that 


CpdP  Sp  dp  0p 

Furthermore,  by  way  of  Eq.  (14)  we  may  write 

50p  =  xSp,  (G6) 

in  which  x  is  evaluated  at  the  original  parcel  position. 

Substituting  Eqs.  (G5)  and  (G6)  into  (G4),  and  then 
Eq.  (G4)  into  (Gl)  leads  to  the  following  expression  for 


the  vertical  buoyancy  acceleration  at  the  new  parcel 
position: 


in  which 


B  =  - 


0p  Sp 


YprclSp, 


SSp/dp  -  X 
Sdp/Sp 


(G7) 


(G8) 


Assuming  that  the  vortex  is  gravitationally  stable,  Yp^;,! 
is  between  zero  and  one;  therefore,  it  generally  reduces 
the  magnitude  of  B.  Y,  defined  in  the  main  text,  is  a 
tempero-azimuthal  mean  of  Yp^j,;. 
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